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Abstract. We present the results of the five mid-IR 
15 /^m (12-18 Aim LW3 band) ISOCAM Guaranteed Time 
Extragalactic Surveys performed in the regions of the 
Lockman Hole and Marano Field. The roughly 1000 
sources detected, 600 of which have a flux above the 80 % 
completeness limit, guarantee a very high statistical sig- 
nificance for the integral and differential source counts 
from 0.1 mJy up to ~5mJy. By adding the ISOCAM sur- 
veys of the HDF-North and South (plus fianking fields) 
and the lensing cluster A2390 at low fiuxes and IRAS at 
high fluxes, we cover four decades in flux from 50 /iJy to 
~0.3 Jy. The slope of the differential counts is very steep 
(a = —3.0) in the flux range 0.4-4 mJy, hence much above 
the Euclidean expectation of a = —2.5. When compared 
with no-evolution models based on IRAS, our counts show 
a factor '^10 excess at 400 /iJy, and a fast convergence, 
with a = —1.6 at lower fluxes. 
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1. Introduction 

Deep galaxy counts as a function of magnitude, or flux 
density, should, in principle, give a constraint on the ge- 
ometry of the universe (llo, Ao). In fact, their departure 
from the Euclidean expectation (no expansion, no curva- 
ture) is dominated by the e-correction (intrinsic evolution 
of the galaxies) and by the /^-correction (redshift depen- 
dence). The understanding of galaxy evolution therefore 
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is a key problem for cosmology, and number counts ap- 
pear to be a strong constraint on the models, which does 
not suffer from the exotic behavior of individual galaxies. 
Most of the energy released by local galaxies is radiated 
in the optical- UV range (Soifer & Neugebauer 1991). If 
this were to remain true over the whole history of the uni- 
verse, then one could follow the comoving star formation 
rate of the universe as a function of redshift by measur- 
ing the optical-UV light radiated by galaxies (Lilly et al 

1996, Madau et al 1996). This scenario changed consid- 
erably after the detection of a substantial diffuse cosmic 
IR background (CIRB) in the 0.1 - 1 mm range from the 
COBE-FIRAS data (Puget et al 1996, Guiderdoni et al 

1997, Hauser et al 1998, Fixsen et al 1998, Lagache et al 
1999) and at 140 - 240 ^^m from the COBE-DIRBE data 
(Hauser et al 1998, Lagache et al 1999). Surprisingly the 
mid-IR/sub-mm extragalactic background light is at least 
as large as that of the UV/optical/near-IR background 
(Dwek et al 1998, Lagache et al 1999), which implies a 
stronger contribution of obscured star formation at red- 
shifts larger than those sampled by IRAS {z > 0.2). To 
understand the exact origin of this diffuse emission and its 
implications for galaxy evolution, we need to identify the 
individual galaxies responsible for it and the best way to 
do that consists of observing directly in the IR/sub-mm 
range. 

In the mid-IR, IRAS has explored the local universe 
(z < 0.2). With a thousand times better sensitivity and 
sixty times better spatial resolution, ISOCAM (Cesarsky 
et al 1996), the camera on-board ISO (Kessler et al 1996), 
provides for the first time the opportunity to perform cos- 
mologically meaningful surveys. Deep surveys have been 
carried out on small fields containing sources well known 
at other wavelengths: the HDF North (Serjeant et al 1997, 
Aussel et al 1999a,b, Desert et al 1999) and the CFRS 
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1452+52 field (Flores et al 1999). This has yielded a small 
but meaningful sample of sources (83 galaxies) with a po- 
sitional accuracy better than 6" , sufficient for most mul- 
tiwavelength studies. Most of these sources can easily be 
identified with bright optical counterparts {Ic < 22.5) 
with a median redshift ofz ~ 0.7 — 0.8 imposed by the 
fc-correction (Aussel et al 1999a, b, Flores et al 1999). Flo- 
res et al (1999) estimate, from their sample of 41 sources, 
that accounting for the IR light from star forming galaxies 
may lead to a global star formation rate which is 2 to 3 
times larger than estimated from UV light only. 

To obtain reliable source count diagrams, better statis- 
tics and a wider range of flux densities are required. 
For this reason, we have performed several cosmological 
surveys with ISOCAM, ranging from large area-shallow 
surveys to small area-ultra deep surveys. These surveys 
were obtained in the two main broad-band filters LW2 
(5-8.5 yum) and LW3 (12-18 /zm), centered at 6.75 /xm and 
15 /um respectively. This paper only presents the source 
counts at 15 /Ltm, because the sample of galaxies detected 
in the 6.75 /^m band is strongly contaminated by Galactic 
stars, whose secure identification requires ground-based 
follow-up observations. Including the surveys over the two 
Hubble deep fields, almost 1000 sources with flirx densities 
ranging from 0.1 mjy to 10 mjy were detected, allowing us 
to establish detailed source count diagrams. Source lists 
for each individual survey, as well as maps and detailed 
description of data reduction, will be found in separate 
papers (see Table |l|) . Forthcoming papers will discuss the 
nature of these galaxies based on the ISOHDF-North sur- 
vey (Aussel et al, in prep. & 1999b), the contribution of 
these galaxies to the cosmic IR background, the relation 
of these observations with ISOPHOT and SCUBA deep 
surveys (Elbaz et al, in prep.), as well as a tentative mod- 
elling of the number counts (Franceschini et al, in prep.). 

2. Description of the Surveys 

The five ISOCAM Guaranteed Time Extragalactic Sur- 
veys (Cesarsky & Elbaz 1996, Table ^ are complemen- 
tary. They were carried out in both the northern (Lock- 
man Hole) and southern (Marano Field) hemispheres, in 
order to be less biased by large-scale structures. These 
two fields were selected for their low zodiacal and cir- 
rus emission and because they had been studied at other 
wavelengths, in particular in the X-ray band, which is an 
indicator of the AGN activity of the galaxies. Only one 
of the 'Marano' maps was scanned at the exact position 
of the original Marano Field (Marano, Zamorani, Zitelli 
1988), while the 'Marano FIRBACK' (MFB) Deep and 
Ultra-Deep surveys were positioned at the site of lowest 
galactic cirrus emission, because they were combined with 
the FIRBACK ISOPHOT survey at 175 /xm (Puget et al 
1999, Dole et al 1999). Indeed the importance of the Galac- 
tic cirrus emission in hampering source detection is much 
larger at 175 /im than at 15 /im, but the quality of the 



two 15 ^m ultra deep surveys in the Marano Field area is 
equivalent. In addition, very deep surveys were taken with 
ISOCAM over the areas of the HDF North (Serjeant et al 
1997) and HDF South. In this paper, we include the HDF 
North results from Aussel et al (1999a), and show for the 
&-st time ISOCAM counts on the HDF South field. 
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A 15, 


fim surveys 
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by increasing 


depth. 












Survey 


Nobs 


Area 830% 


tint 


#gal 


Slope 


Name 




(am ) (mjy) (min) 






(1) 


(2) 


(3) (4) 


(5) 


(6) 


(7) 


Lockman Shallow^'"' 


3 


1944 1 


3 


80 


-2.1 ±0.2 


Lockman Deep'"' 


6 


510 0.6 


11 


70 


-2.4 ±0.3 


MFB Deep(") 


18 


710 0.4 


15.4 


144 


-2.4 ±0.2 


Marano UD^"' 


75 


70 0.2 


114 


82 


-1.5 ±0.1 


MFB UD'") 


75 


90 0.2 


114 


100 


-1.5 ±0.2 


HDF North*'') 


64 


27 0.1 


135 


44 


-1.6 ±0.2 


HDF South*") 


64 


28 0.1 


168 


63 


-1.4 ±0.1 


A2390''^) 


100 


5.3 0.05 


432 


31 


-1.2 ±0.6 



Comments: Col.(l) Survey name with reference: (a) in prepa- 
ration, (b) Aussel et al (1999), (c) Altieri et al (1999); Col.(2) 
maximum number of pointings on the same sky position (re- 
dundancy); Col. (3) the total area covered in square arcminutes; 
Col. (4) the flux at which the survey is at least 80% complete; 
Col. (5) the corresponding integration time per sky position (in 
minutes); Col. (6) the number of galaxies whose flux is over the 
80% completeness threshold; Col. (7) the slope of the fit to the 
integral log A'' — log S. A2390 completeness limit includes the 
corrections for lensing magnification. 



3. Data Reduction and Simulations 

The transient behavior of the cosmic ray induced glitches, 
which makes some of them mimic real sources, is the 
main limitation of ISOCAM surveys. We have developed 
two pipelines for the analysis of ISOCAM surveys in or- 
der to obtain two independent source lists per survey 
and improve the quality of the analysis. PRETI (Pattern 
REcognition Technique for ISOCAM data), developed by 
Starck et al (1999), is able to find and remove glitches 
using multi-resolution wavelet transforms. It includes also 
Monte Carlo simulations to quantify the false detection 
rate, to calibrate the photometry and to estimate the com- 
pleteness. The 'Triple Beam-Switch' (TBS) technique, de- 
veloped by Desert et al (1999), treats micro-scanning or 
mosaic images as if they resulted from beam-switching ob- 
servations. All the surveys have been independently ana- 
lyzed using both techniques and the source lists were cross- 
checked to attribute quality coefficients to the sources. 
PRETI and TBS agree at the 20 % level in photometry 
(corresponding to the photometric accuracy of both tech- 
niques), and with an astrometric accuracy smaller than 
the pixel size (due to the redundancy). PRETI allowed 
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Fig. 1. Integral counts, i.e. the number of galaxies, N, de- 
tected at 15 /im above the flux S(mJy), with 68 % con- 
fidence contours. K counts (Gardner et al 1993) and B 
counts (Metcalfe et al 1995), multiplied by the ratio v jvx^ 
to represent the relative energy densities at high fluxes, 
are overplotted with open circles and filled squares, re- 
spectively. The hatched area materializes the range of pos- 
sible expectations from models assuming no evolution and 
normalized to the IRAS 12 /im local luminosity function 
(LLF). The upper limit was calculated on the basis of 
the LLF of Rush, Malkan & Spinogho (1993), as Xu et 
al (1998) and shifted from 12 to 15 /im with the template 
SED of M82; the lower hmit uses the LLF of Fang et al 
(1998) and the template SED of M51. 
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Fig. 2. Differential Number Counts of 15 /im Galaxies, 
with 68% error bars. The counts are normalized to a Eu- 
clidean distribution of non-evolving sources, which would 
have a slope of a = —2.5 in such a universe. Data 
points: A2390 (open stars), ISOHDF-North (open circles), 
ISOHDF-South (filled circles), Marano FIRBACK (MFB) 
Ultra-Deep (open squares), Marano Ultra-Deep (crosses), 
MFB Deep (stars), Lockman Deep (open triangles), Lock- 
man Shallow (filled triangles). The hatched area materi- 
alizes the range of possible expectations from models as- 
suming no evolution (see Fig.l). 



us to attain fainter levels in deep surveys, whereas in the 
shallow surveys, where the redundancy is not very high, a 
very strict criterion of multiple detections had to be ap- 
plied. Finally, we have made Monte Carlo simulations by 
taking into account the completeness and the photometric 
accuracy to correct for the Eddington bias and to compute 
error bars in the number count plots. 

4. The ISOCAM 15 /xm source counts 

Figures |l| and g show respectively the integral and the 
differential number counts obtained in the five indepen- 
dent guaranteed time surveys conducted in the ISOCAM 
15 /xm band, as well as the HDF surveys. The contribution 
of stars to the 15 /xm counts was corrected. It is negli- 
gible at fluxes below the mJy level as confirmed by the 
spectro-photometric identifications in the ISOHDF-North 
(1 star out of 44 sources). South (3 stars over 71 sources), 
and CFRS 1415-1-52 (1 star over 41 sources ranging from 
~ 0.3 mJy to ~ 0.8 mJy). In the Lockman Shallow Sur- 
vey {S\^p,rn > 1 mJy), about 12 % of the sources were 
classified as stars from their optical-mid IR colors (us- 
ing the Rayleigh- Jeans law) . We have also represented the 
counts from the ISOHDF-North (from Aussel et al 1999a), 
ISOHDF-South, and, at the lowest fiuxes, the counts ob- 



tained from the A2390 cluster lens (down to 50/iJy, in- 
cluding the correction for lensing magnification; Altieri et 
al 1999, see also Metcalfe et al 1999). We have only in- 
cluded the flux bins where the surveys are at least 80 % 
complete, according to the simulations. 

The first striking result of these complementary source 
counts is the consistency of the eight 15 /im surveys over 
the full flux range. Some scatter is nevertheless apparent; 
given the small size of the fields surveyed, we attribute it 
to clustering effects. The two main features of the observed 
counts are a significantly super-euclidean slope (a = —3.0) 
from 3 to 0.4 mJy and a fast convergence at flux densities 
fainter than 0.4 mJy. In particular, the combination of flve 
independent surveys in the flux range 90-400 /iJy shows 
a turnover of the normalized differential counts around 
400 /iJy and a decrease by a factor ~ 3 at 100 /iJy. We 
believe that this decrease, or the flattening of the integral 
counts (see the change of slope in col(7) of Table ^ below 
~400 /iJy, is real. It cannot be due to incompleteness, since 
this has been estimated from the Monte-Carlo simulations 
(see Section ||) . The differential counts can be fitted by two 
power laws by splitting the flux axis in two regions around 
0.4 mJy. In units of mjy~^ deg~^, we obtain, by taking 
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into account the error bars (S* is in mJy): 

(2000 ±600) S'(-i'^±°-2) ... 0.1<S'<0.4 



dNjS) 
dS 



(470 ±30) S'^-^oioi) ... 0.4<5<4 



(1) 



In the integral plot, the curves are plotted with 68 % confi- 
dence contours based on our simulation analysis. The total 
number density of sources detected by ISOCAM at 15 /im 
down to 100 /iJy (no lensing) is (2.4 ± 0.4) arcmin"^. It 
extends up to (5.6 ± 1.5) arcmin"^, down to 50/iJy, when 
including the Icnsed field of A2390 (Ahieri et al 1999). 

5. Discussion & Conclusions 

We have presented the 15 /im differential and integral 
counts drawn by several complementary ISOCAM deep 
surveys, with a significant statistical sampling (993 galax- 
ies, 614 of which have a flux above the 80 % completeness 
limit) over two decades in flux from 50 /iJy up to 5 mJy. 

The differential counts (Fig. |2[), which are normalized 
to S^'^'^ (the expected differential counts in a non expand- 
ing Euclidean universe with sources that shine with con- 
stant luminosity), present a turnover around 5*15 ^m=0. 4 
mJy, above which the slope is very steep {a = — 3.0±0.1). 
No evolution predictions were derived assuming a pure 
fe-correction in a flat universe (go = 0-5), including the ef- 
fect of Unidentified Infrared Bands emission in the galaxy 
spectra. In the Fig.|l|&^, the lower curve is based on the 
Fang et al (1998) IRAS 12/j,m local luminosity function 
(LLF), using the spectral template of a quiescent spiral 
galaxy (M51). The upper curve is based on the Rush, 
Malkan & Spinogho (1993) IRAS-12/xm LLF, translated 
to 15 /im using as template the spectrum of M82, which is 
also typical of most starburst galaxies in this band. More 
active and extincted galaxies, like Arp220, would lead to 
even lower number counts at low fluxes while flatter spec- 
tra like those of AGNs are less flat than M51. In the ab- 
sence of a well established LLF at 15 fim, we consider these 
two models as upper and lower bounds to the actual no- 
evolution expectations; note that the corresponding slope 
is ~ —2. The actual number counts are well above these 
predictions; in the 0.3 mJy to 0.6 mJy range, the excess is 
around a factor 10: clearly, strong evolution is required to 
explain this result (note the analogy with the radio source 
counts, Windhorst et al 1993). 

We believe, according to the results obtained on the 
HDF and CFRS fields (Aussel et al 1999a,b, Flores et al 
1999), that the sources responsible for the 'bump' in the 
15 /im counts are not the faint blue galaxies which domi- 
nate optical counts at z ~ 0.7. Instead, they most proba- 
bly are bright and massive galaxies whose emission occurs 
essentially in the IR and could account for a considerable 
part of the star formation in the universe at z < 1.5. 

In Fig. ^ we have overplotted the integral counts in 
the K (Gardner et al 1993) and B (Metcalfe et al 1995) 
bands, in terms of i^Si,. For bright sources, with densities 



lower than 10 deg~^, these curves run parallel to an in- 
terpolation between the ISOCAM counts presented here 
and the IRAS counts; the bright K sources emit about 
ten times more energy in this band than a comparable 
number of bright ISOCAM sources at 15 /tm . But the 
ISOCAM integral counts present a rapid change of slope 
around 1-2 mJy, and their numbers rise much faster than 
those of the K and B sources. The sources detected by 
ISOCAM are a sub-class of the K and B sources which 
harbor activity hidden by dust. Linking luminosity to dis- 
tance, we predict a rapid change of the luminosity function 
with increasing redshift, which can only be confirmed by 
a complete ground-based spectro-photomctric follow-up. 
We should be able to follow the evolution of the luminos- 
ity function from z ~ 0.2 to 1.5 with the large number of 
galaxies detected in the Marano Field surveys. The combi- 
nation of the intensity of the Ha emission line (redshifted 
in the J-band) with the IR luminosity could set strong 
constraints on the star formation rate. Finally, emission 
line diagnostics, combined with hard X-ray observations 
with XMM and Chandra, would allow us to understand 
whether the dominant source of energy is star formation 
or AGN activity. 
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